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■ Abstract. A direct procedure for determining the propagator associated with a quantum 

mechanical problem was given by the Path Integration Procedure of Feynman. The Green 
• function, which is the Fourier Transform with respect to the time variable of the propagator, 

can be derived later. In our approach, with the help of a Laplace transform, a direct way to 
get the energy dependent Green function is presented, and the propagator can be obtained 
later with an inverse Laplace transform. The method is illustrated through simple one 
dimensional examples and for time independent potentials, though it can be generalized to 
O" 1 the derivation of more complicated propagators. 
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'■ 1 Introduction 



It is well know that quantum mechanics acquired its final formulation in 1925-1926 through 
fundamental papers of Schrodinger and Heisenberg. Originally these papers appeared as two 
independent views of the structure of quantum mechanics, but in 1927 Schrodinger established 
their equivalence, and since then one or the other of the papers mentioned have been used to 
analyze quantum mechanical systems, depending on which method gave the most convenient 
way of solving the problem. Thus the existence of alternative procedures to solve a given problem 
can be quite fruitful in deriving solutions of it. 

In the 1940's Richard Feynman, and later many others, derived a propagator for quantum 
mechanical problems through a path integration procedure. In contrast with the Hamiltonian 
emphasis in the original formulation of quantum mechanics, Feynmans approach could be re- 
ferred to as Lagrangian and it emphasized the propagator K(x,t,x' ,t') which takes the wave 
function tp(x', t') at the point x' and time t' to the point x, at time t, i.e. 

i>{x,t)= I K{x,t,x',t')-ifj(x',t')dx'. (1) 



While this propagator could be derived by the standard methods of quantum mechanics, 
Feynman invented a procedure by summing all time dependent paths connecting points x, x 1 and 
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this became an alternative formulation of quantum mechanics whose results coincided with the 
older version when all of them where applicable, but also became relevant for problems that the 
original methods could not solve. Feynmans procedure first led to the propagator K(x,t,x' ,t') 
and then by a Laplace transform to the corresponding Green function G(x, x' , E) with E being 
the energy. We found Feynmans method for deriving the propagator, though entirely correct, 
somewhat cumbersome to use, and thus tried to look for alternative procedures. 

As we mentioned before, in Feynmans approach the first step is deriving the propagator 
K(x,t, x' ,t') and later the energy dependent Green functions G(x,x' , E). In this paper we 
invert the procedure, we start by deriving the G(x,x' , E) which is a simpler problem, at least 
in the one dimensional single particle case we will be discussing here. Once we have G(x, x' , E) 
the K(x,t,x' ,t') is given by the inverse Laplace transform and can be written as 



where c is a constant that allows the upper line ihc + E in the complex plane of E to be above 
all the poles of G(x, x' , E). For compactness in the notation from now on we will take t' = so 
we write K(x,t,x' ,t') as K(x,x',t). 

The real hard part in our approach will be the determination by ([2]) of K(x, x' , t) but this is 
a well defined problem in mathematics and procedures have been developed to solve it. 

This is then the program we plan to follow. In Section [2] we show that for a single particle 
in one dimension (the initial case of our analysis) all we need to know are two independent 
solutions itg of the equation. 



to be able to derive G(x, x' , E) in Section [3j 

We then consider in Section HJ three elementary cases, the free one dimensional particle, the 
corresponding one with a 5 interaction at the origin x = 0, and the harmonic oscillator. In the 
first two cases the integral ([2]) is trivial to evaluate. In the case of the harmonic oscillator the 
evaluation of ([2]) requires a more careful analysis but it can be carried out. In all three cases 
our final result is identical to the one presented in the book of Grosche and Steiner PQ that use 
Feynmans method to derive the results. 

Thus we have an alternative method for deriving K (x, x', t) but it remains to be shown that 
it can be applied to more particles in more dimensions and of arbitrary angular momenta and 
whether the analysis can be extended to relativistic as well as time dependent problems. 

What ever results may be obtained in the future, it seems that our alternative approach 
follows more closely the standard procedures of quantum mechanics and could be useful in 
a simpler derivation of the propagators. 

2 The Hamiltonian of the problem and the equation 
for the propagator 

We start with the simplest Hamiltonian of one particle in one dimension, i.e. 




(2) 



\-h 2 d 2 





with thus far an arbitrary potential V(x). 
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From the equation dU) that defines the properties of the propagator it must satisfy the 
equation 

h 2 d 2 d~ 

K(x,x',t) = (3) 



2m dx 2 ^ ^ dt 



and besides if t = it becomes 
K(x, x , 0) = 5(x — x). 
We proceed now to take the Laplace transform of (J3j) 



f 

Jo 



exp(— st) 



h2 92 „, \ -, d ' 

2^d^ 2+v{x) - lh m 



K(x,x',t)dt 



h 2 d 2 G{x,x',s) 
2m dx 2 



+ V(x)G(x, x', s) - ih r exp{-st) dK ^ X ' X,,t ^ dt = 0, 
Jo ot 



where 



G(x,x',s) = / e- st K(x,x',t)dt. 
Jo 

We note though that 



f 

Jo 



exp(-st) dK ^ X ' ,t ^ dt= ^-[e- st K(x,x',t)]dt + s f°° e - st K(x,x',t)dt 
= —5{x — x') + sG(x, x' , s), 



where we made use of and ([5]). 

With the help of (|6|) we see that G(x,x',s) satisfies 



2m dx 2 



+ V(x) — ihs 



G(x, x', s) = —ih5{x — x'), 



(4) 



(5) 



(6) 



(7) 



where we now have that the partial derivative with respect to x becomes the ordinary one as 
there is no longer a time variable. 

We integrate ([7]) with respect to x in the interval x' — e < x < x' + e and in the limit e — > 
obtain two equations 



ff_ /dG 

2m \ dx 

h 2 d 2 



h 2 (dG 

x =x'+o V. dx j x=x ,_q 



-ih, 



2m dx 2 



+ V(x) — ihs 



G(x,x',s) = 0, x^x'. 



(8) 
(9) 



We proceed now to indicate how we can derive the explicit expression of K(x, x', t) with the 
help of the Green function G(x,x',s) of the corresponding problem satisfying (jSJ) and Q. 



3 Determination of the Green function 

and the inverse Laplace transform for the propagator 

Our interest is not to stop at equations ([8]), ([9]) for G(x,x',s) but actually to get K(x,x',t) for 
which we can use the inverse Laplace transform [2] to get 



K(x, x', t) = / G(x, x', s)e st ds, 

2vri ./, 



(10) 
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where the integration takes place along a line in the complex plane s parallel to the imaginary 
axis and at a distance c to it so that all singularities of G(x, x' , s) in the s plane are on the left 
of it. 

To have a more transparent notation rather than the s plane we shall consider an energy 
variable E proportional to it through the relation 



or 



-i(E/h) 



E = ihs 
and define G(x,x',E) by 

-%G{x, x', E) = G(x, x', -iE/h). 
The energy Green function must be symmetric under interchange of x and x', i.e. 
G{x,x',E) = G{x',x,E) 
which combines with the two equations ((8|), (|9|) to give in this notation 



dG~ 




'dG~ 




2m 


dx 


x=x'+0 


dx 


x=x' — 


h 



2m dx 2 



+ V{x) - E 



G(x,x',E) = for x^x'. 



(11) 

(12) 
(13) 



Let us first consider the case when x < x' and proceed to show that the equations (Ilip - (ll3p 
determine in a unique way the Green function of the problem. For this purpose we introduce 
with the notation u E (x) two linearly independent solutions of the equation f) 1 31) 



h 2 d 2 



+ V(x) - E 



2m dx 2 

From this equation we see that 



u%{x) = 0. 



u E {x 



d 2 u~t(x) ,d 2 Up(x) d 

~ ut = -r 

ax dx dx 



- du E 
' dx 



ul 



dx 



0. 



Thus the Wronskian of the problem defined by 
W{E)=u E {x)-^-u E {x)- 



dx 



(14) 



is independent of x. 

As G(x,x' ,E) satisfies (|13h we can write it for x < x' as 

G(x,x',E) = F(x',E)u E {x), 



(15) 



choosing one of the two solutions of equation (|13h and F(x', E) is as yet an undetermined 
function of x', E. 

We see from the symmetry of G(x,x' , E) that it must satisfy the same equation (|13j) in x' so 
that from (|15[) we get 



h 2 d 2 



+ V{x) - E 



F(x',E) = 



2m dx' 2 

and thus F(x',E) must a be linear combination of the two independent solutions u E (x), i.e. 
F(x',E) = a + {E)u E (x') + a-(E)u E (x') 
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and our Green function becomes 

G(x,x',E) = [a + {E)u + E (x') + a.(E)u-{x')]u + E (x), (16) 
while for the other case, i.e. x > x', the symmetry of the Green function demands 

G(x, x', E) = [a + (E)u E {x) + a_(E)u E {x)] u%{x'). (17) 
Replacing (fT6j) and (fT7|) in (fT2|) we find that the coefficient a + {E) vanishes and a-{E) satisfies 

2 777 

a-(E)W(M) = --j-. (18) 

Thus from <$ja§, dTTJ) and CEH]) we get that 

, , 2m i , . I u7?{x')ut(x) if x < x', 
G(x,x',E) = -—W- 1 (E)\ I)'?;' ' (19) 

n u E {x)u E {x ) it x > x . 

We thus have the explicit Green function of our problem once we can obtain two independent 
solutions of the equations (fT3j) . 

Once G(x,x',E) has been determined, the propagator K(x,x',t) is given by the inverse 
Laplace transform (|10p which in terms of the E variable becomes 

i rihc+oo 

K(x,x',t) = — — / exp(-iEt/h)G(x,x',E)dE, (20) 
2vrfe J^.oo 

where now the integral takes place in the E plane over a line parallel to the real axis with all 
the poles of G(x,x',E) below it. 

We proceed to give some specific examples of application of our method. 

4 Specific examples 
a) The free particle 

The potential V{x) is taken as zero and so the equation (fT3|) becomes 

h 2 d 2 



E 



G{x,x',E) = 0. 



2m dx 2 

We introduce the variable k through the definition 

^ h 2 k 2 lrn h 2 k „ 

E = , dE = dk 

2m m 

and thus the Ug(x) for this problem satisfy the equation 

~d 2 l2 
+ k 2 



dx 



u E {x) = 0, u E (x) = exp(±ikx) 

with the Wronskian (I14p given by 
W(E) = 2ik. 

Thus from the two cases of (|19|) our function G(x, x' , E) is written compactly as 

G(x,x'E) = — exp[ik\x-x'\}. (21) 
nk 
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The propagator K(x,x',t) is given by (|20|) in terms of G(x,x',E) and substituting (|2"T|) in it 
and writing in terms of k we get 

1 1'°° 

K (x, x , t) = — / exp[ik\x — x' | - i(hk 2 /2m)t]dk, (22) 

where, as G(x,x' , E) has no singularities, the energy can be integrated over the real line 
— oo < E < oo while k has double the range of E. 

The integral (|22p can be determined by completing the square and we get 

im(x-x') 2 1 ^ 



K(x,x',t) = \/^~^J ex P 



2ht 



which has been derived also by many other methods, 
b) The case of the delta potential 

We wish now to discuss the effect on the Feynman propagator of a potential 
V(x) = Q(x) + b5(x), 

where Q{x) in a continuous function of x and we assume b > to avoid bound states of the 5 
potential. 

The equation for Ug becomes now 

[H + bS(x) -E]u%{x) = 0, 

where 

h 2 d 2 

2m dx 2 ^ 

The Green function G(x, x' , E) satisfies the equations (|12|) . (|13|) which can be written as the 
single equation 

[H + bS(x) - E]G{x, x', E) = -iH5(x - x') (24) 

and (|12p holds if we integrate (|24p with respect to the variable x in the integral x' — e<x<x' + e 
in the limit e — > 0, and the one corresponding to (113j) holds when x ^ x' . 

For i/Owe have (|24|) with no delta potential and therefore G can be written as 

G(x,x',E) = G Q (x,x',E) +F(x,x',E), x ^ 0, (25) 
where Gq(x,x' , E) is the Green function satisfying 

[H -E]Gq{x,x',E) = -ih5{x-x) (26) 
while F(x,x' , E) is a solution of the corresponding homogeneous equation, i.e. 

[H - E]F(x, x', E) = 0, x^x'. 

and the form of F(x,x' ,E) is to be determined. The continuity of G(x, x' , E) at x = allows 
to write (|25p for all values of x, x' and with this in mind we can replace (|25l) in (|24h to obtain 

[H-E- b5(x)} [G Q (z, x, £?) + F(x, F)] = -iM(x - x') 

= -iM(x - x') - 65(x)G Q (0, x', E) + [H-E- b5{x)]F{x, x' , E), (27) 
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where in the second line we have used (|26|) . The two lines in f)2T|) imply 

[H-E- bd(x)]F(x, x',E) = bS(x)G Q (0, x' , E) 
which is a version of (|24p but with a source term. Therefore the solution can be readily given as 

■ poo 

F(x,x',F) = — — / dx"G(x,x",E)(bd(x")G Q {0,x',E)). 

" J —oo 

Performing the integral in the last expression and using f)25[) for G, we have 

F(x, x r , F) = 7 [G Q (x, 0, F) + F(x, 0, E)] G Q (0, x', E), (28) 

where 7 = —ib/h. To determine F(x,0,F) in the RHS of ([28]) we set x' = and solve for F, 
obtaining 

_ 7 G Q (x,0,F)G Q (0,0,F) 

(x,u ' j " i- 7 Gq(o,o,f) ' l29j 

Finally, (I29p can be replaced back in (|28[) to get 

, m 1 G Q (x,0,E)G Q (0,x',E) 
F ^ X ^= 1- 7 G Q (0,0,F) • (30) 

With this, G(x, x', F) is given now in terms of Gq(x, x' , E) and if we make Q = we can apply 
it to the case of the free particle. The Green function Gq(x, x' , E) is given by (|21|) and the Green 
function of our problem becomes 

, im r .,, , n m 2 bexp[ik(\x\ + \x'\)} 
G(x, x , F) = — exp[^|x - x |] - ^ , „ , ;mb . 



hk" ~ IJ 2fr4 + ^ 

im . , . ... exp[i/c(|x| + Ix'l)] im exp[i/c(|x| + \x'\ 
■ exp[ik\x - x'\] - —2 1 — + 



"nib L 1 IJ 2fr 2 it 2ft 2 A; + i# 

This is the same result that appears in Grosche and Steiner [TJ (6.12.4), p. 328] and accounts ex- 
plicitly for the four possible cases ±x, ±x' > 0. The inverse Laplace transform of this expression 
can be easily evaluated since it contains integrals of the form 



/oo 
dkk~ n exp (-Ck 2 + Dk) , 
-00 



where C and D are independent of k and the integrals are either gaussians or error functions 
when n = or n = 1 respectively. The propagator is then 

mb ( mb ( ' , ,, imb 2 t 



1 1 m b ( m b ( \ i\ imb t\\ 

K(x, x , t) = K (x, x , t) + ^2 ex P ( ~~fi2 I M + F I + -^3- J J 



, m / ,,,,/, i&i 
x erfc < \ / — — x + x 



2ifa V 1 1 1 ft 

with Kq(x,x' ,t) as in (|23p . This result coincides with the one that Grosche and Steiner [TJ 
(6.12.2), p. 327] obtained by the Path Integral Methods. 

Note also that if in G of (|25|) we replace the F(x,x',F) by (|3U|) we also get the result of 
Grosche and Steiner [TJ (6.12.1), p. 327]. 
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c) The harmonic oscillator 

The potential V(x) is proportional to x 2 and thus Ug(x) satisfies the equation 

h 2 d 2 1 



2m dx 2 2 



+ -muj 2 x 2 — E 



u E (x) = 0, 



(31) 



where uj is the frequency of the oscillator. 
We introduce the variables 



2muj 



n 



E 1 
V ~hL~2 



in terms of which the equation ()3ip takes the form 



d 2 z 2 1 



u%{x) = 0. 



(32) 



(33) 



Two independent solutions of (|33p are given by parabolic cylinder functions [3], i.e. 
u%{x) = D p (±z). 

The Wronskians of these functions, where the derivative is taken with respect to the x rather 
than the z variable, is from (|14p and (|32p given by 



W{E) 



2moj 



2mw 



Dp(-z) 



D p (-z) 



dD p (z) 
dz 



D p {z) 



dD p {-z) 



dz 



1 



■D p+ x(z) + -zD p (z) 



+ D p (z) 



D p+1 (-z) - -zD p (-z) 



h 

where we made use [U (9.247-3), p. 1066] together with the fact that 

[dD p (-z)/dz] = -[dDp{-z)/d{-z)\. 

The Wronskian is independent on z so we may take any value of the latter and we choose 
z = to get 



W(E) = -^^2Dp(0)D p+1 (0). 

We note from [31 (9.240), p. 1064] that we can write D p (z) in terms of the degenerate 
hypergeometric function $ and, in particular 



D p (0) = 2^-^-$ (*Lo) 
PK ' T(k^) V 2' 2' J 



while from [3 (9.210), p. 1058] the $(-§, 4,0) = 1 so that finally 



W(E) 



2muo 2' p+1 tt 

~r(i?)r(-f) 



2muj \/tt 

h v(- P y 



(34) 



where for the last expression in (|34|) we made use of the doubling formula for the T function 
given in O (8.335), p. 938]. 
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From the general relation (|19p we then obtain that the Green function of the oscillator is 
given by 



G(x,x',E) = ^—r(-p)D p (z)D p (-z') (35) 

with p, z given by (|32p . z < z' and z 1 has the same definition as z but x replaced by x' . When 
we consider the case z > z' we have an expression which is similar to ([35|) but interchanging z 
and z' . The complete formula can be written in a compact way introducing the variables 

x> = max{x,x'}, x < = min{x,rr'} 

and thus we have 



G(x,x',E) = J^-T ^ - i-) Dei J*£L x> D,j ( -d^ X< ) . (36) 



irhw \2 hw J hui 2 \ V h 



hu 2 




The expression (j36j) is identical to [H (6.2.37), p. 179] except for a factor of \/2. 

We want though to obtain K(x,x',t) of (|20p . but for this we can analyze the pole structure 
of (1361) in order to evaluate the inverse Laplace transform of G(x, x' , E) by means of the residue 
theorem. The result of this procedure is the series known as the spectral decomposition of 
K(x,x',t). This series can be evaluated by using an identity of Hermite polynomials in [3] and 
references cited therein. We carry out the analysis in the Appendix and get the final result 



. , . / muj \ 1 / 2 f imoj 
Kix > X > t) = \ 2*ihsmJ eXP \ 



2hs'm ujt 

which coincides with expression in [U p. 160]. 



(x + x ) cos cut — 2xx' 



5 Conclusion 

We will state here the full steps to get the Feynman propagator of a non-relativistic single 
particle one dimensional problem. 
Our Hamiltonian is 

h 2 d 2 



H 



2m dx 2 



+ V(x) 



with an arbitrary potential V(x). 

We assume that two independent eigen-functions of energy E can be found and denoted 
by u%{x). 

We need then to determine the Wronskian 

du~^(x) +f ^du~(x) 



W(E)=u-(x)—^-u+(x) 



dx dx 

whose value, as we know, is independent on x. 

Following the analysis of Section [3] we can determine the energy Green function 

G(x, x',E) = -^W{EY l { u e( x >e( x ) iix < x 'i 
h I u E (x)ug(x') if x > x'. 

Using the Laplace transform the Feynman propagator becomes 

i rihc+oo 

K(x,x',t) = — -/ exp(-iEt/h)G(x,x',E)dE. (37) 
2nhi J ih c-oo 



10 



M. Moshinsky, E. Sadurm and A. del Campo 




Figure 1. Integration contour for the determination of the propagator. 

Once we get u E (x) the only troublesome part of our calculation is the integral (|37p as we 
already saw in the discussion of the examples in Section |H 

For time dependent Hamiltonians probably other techniques should be used but we have not 
developed them yet. 

A Determination of the propagator for the harmonic oscillator 

We shall proceed to express the Green function G(x,x' , E) in terms of its energy poles and 
residues and then apply the Laplace transform to it to get the propagator. 
To achieve our objective we start with the complex variable integral 

I e ^{pH>AT(-p)D p (z)D p (-z)dP, (38) 
J c 

where 



E 1 2muo ! 2moj , / 2m 

p= j^-r z = y— x > z = y— x > A= y^ (39) 

and the contour C of integration is depicted in Fig. [H 

The function T is gamma while the parabolic cylinder is function D p (z). If we substitute p, 
z, z' of (f39j) in the integration (i38j) and replace the contour in Fig. [Qjust by its upper line we 
obtain the integral 

1 rihc+co j 2 m / 1 e \ 

K(x,x',t') = — -/ exp(-iEt/h)\ /-= — r[ -- — ]£>« i (z)D e i(-z')dE. 

We proceed now to evaluate integral (|38p in terms of energy poles and residues of its integrand. 
The parabolic cylinder functions D p (z), D p (—z') are analytic in the full p plane as shown in 
[31 (2), p. 1066], while T(—p) has poles only on the real p axis with values 

p = n = 0,1,2,.... (40) 

Note also that the residues of r(— p) at the poles given in (f40l) are 

Resr(-p) = ^-^ 
n\ 

as shown at |3 p. 933]. 
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Furthermore the integral 
1 f°° D p {z)D p {-z>) 



2m 



J-oo P~ n 



as D p (z) is analytic in the full p plane. 

The lower circle in the contour of Fig. Q] does not contribute to the integral because the term 
exp(— iutp) vanishes when p — > oo with a negative imaginary part. As the poles marked with 
across in Fig. Q] are only on the real p axis we can eliminate the region between the upper line 
and the real axis and all that remains is the residue of T(—p) in the integrand (I38p so K(x, x' , t) 
takes the form 



The parabolic cylinder of index n can be put in terms of an Hermite polynomial 

D n {z) = 2-5 e"T H n (-?L 
so the propagator can also be written as 



, . / 2mw \ - e 



hir ^ n\2 n 

n=0 



xexp{-^(, 2 + x' 2 )}F n (^x)^ 

Using the relation [5] 

™ H n (z)H n (z') 1 _j2^(z 2 + z' 2 )-2zz') 

2-^ n\ ** 



IT* 



7r=1F exp l « 

we obtain 



n 

n=0 



K(x,x';t) = Jj—e ~ exp [- — (x 2 + x' 2 )j 
x exp 



V 7 ! - e~ i2ujt 
muj e- iujt {e- iujt (x 2 + x' 2 ) - 2xx') 



2mw 1 [ mu> 

exp 



[X + X 



e~ 2iujt 1\ 2e~ iujt xx' 



irh yje iu)t — e~ iujt [ h 

I muj f imuj r . , . 9 , 9 . n ] . , , 

= J exp<^— ; — - cos (cot )(x 2 + x /2 )- 2xx'} }. (41) 

As equation (|4ip is identical to (|37p and this also applies to the other corresponding equations 
discussed in the examples, we see that our alternative approach is equivalent to that of Feynman. 

Other methods for deriving the propagator for the case of the harmonic oscillator have been 
proposed recently [51 E] . 
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